General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



N83-20952 


(NASA-CE“ 170118) EnHAn C eD HANUAL 
CONTROLLABILITY VIA ACTIVE CONTROL OR 
AEROELASTIC VEHICLES Semiannual Status 
Report, 1 Aug. 1982 - 31 Jan. 1983 (Purdue 
Univ.) 58 p liC A04/MF A0 1 CSCL 01C G3/08 


Unclas 

03226 


ENHANCED MANUAL CONTROLLABILITY 
VIA ACTIVE CONTROL OF 
AEROELASTIC VEHICLES 


1 


& 










iP 


4 


Semi-Annual Status Report 
For the Period Aug. 1, 1982 to Jan. 31, 1983 


Grant No. NAG-1 -254 


Principal Investigator: Dr. David K, Schmidt 

School of Aeronautics & Astronautics 
Purdue University 
W. Lafayette, IN 47907 


NASA Technical Officer: Mr. William Grantham 

Theoretical Mechanics Branch 
M.S. 152A 

NASA Langley Research Center 
Hampton, VA 23665 


March 30, 198-T 




i 

i 

f 

i 

\ 

| r 

l 

i 

i 

i 


Page 

1 

1 

22 

24 

40 


Content s 

1 . Introduction 

2. Technical Progress 

3. Personnel and Plans 
Appendix A Extended Abstract 
Appendix B Longitudinal Residues 


1. Introduction 


This constitutes the semi-annual progress report for the period 
7/31/82-1/31/83 on the research being performed by the School of Aeronau- 
tics and Astronautics, Purdue University, for NASA's Langley Research 
Center under grant number NAG-1 -254. The subject area of the research 
includes flexible aircraft dynamics, handling qualities of flexible air- 
craft, and active control for good handling characteristics for such 
vehicles. 

2. Technical Progress 

During this first year of grant activity, a modal analysis technique 
has been developed for evaluating the effects of elastic modes on aircraft 
dynamic response, and the handling qualities implication of these effects. 
This technique will be discussed, along with results of application of the 
technique in a technical paper entitled 

A Modal Analysis of Flexible Aircraft Dynamics With 
Handling Qualities Implications . 

This paper will be presented at the 1983 AIAA Atmospheric Flight Mechanics 
Conference, to be held in Gatl inburg, Tenn. in August. The extended 
abstract for this paper is included as Appendix A of this status report. 

As discussed in the paper's abstract, the analytical technique for 
the analysis of flexible flight vehicle dynamics includes the analysis 
of the vehicle's dynamic modes, and specifically includes the dynamic 
model of the atmospheric disturbance and pilot lags in the analysis. 
Fundamental results obtained in such an analysis are not only the eigen- 


values and eigenvectors for the system modes, but also the degree of 
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controllability (disturbabil ity) , o bservability , and impulse response 
residues for each system mode. By obtaining these important measures of 
response to meaningful system inputs, such as those from the pilot and 
turbulence, the significance of the participation of each mode may be 
assessed. The analytical methodology is briefly summarized in Table I 
included here. 

As of the date of this report, the longitudinal dynamics of a family 
of eight vehicle math models have been evaluated, and the lateral direc- 
tional dynamics are under investigation presently. These eight configu- 
rations were evaluated by Yen in a fixed-base laboratory simulation study. 
The geometry of all of these configurations is similar to the B-l air- 
craft, but the structural stiffness of each configuration is varied para- 
metrically in such a way that the elastic mode shapes remain unchanged, 
but the in-vacuo frequencies are reduced. These eight configurations 
then are identified in terms of their eigenvalues in Table II. Shown in 
Figs. 1 and 2 are the attitude tracking results and subjective ratings 
obtained in Yen's simulations. 

Note, for example, that in terms of eigenvalues alone, it is not 
at all clear why the results for configuration 3 were so much worse than 
for configuration 4. Both had unstable phugoid roots, and one rather low- 
frequency aeroelastic mode. In the case of configuration 3, this low- 
frequency mode included a large amount of fuselage bending, while for 
configuration 4, the eigenvectors indicated that the low-frequency mode 
included a good deal of wing bending. This fact partially explains why 
attitude tracking would be more difficult for configuration 3. 
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SUMMARY OF MATH METHOD 
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and outputs (responses) 


y = 0x v + Ex p + Fx g 


= [E F C ] 


Table I (cont'd) 


so collecting we have 




X a A*X + G*to + G*to„ 
p p g 9 

y = C*x 


Now Define 

T ■ modal matrix (matrix of eigenvectors) of A* 
Then in Terms of Vehicle Modal States 
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B = [b 1 |b 2 {•••■ b^,] ; D * [d] 


Note for each output y 

y j (s> . ? 

v 75 ” \ 




each input we would have 
; = Residue for 


So the Residues Are For Each y. and u or to are 

1 p j 9 


R ij “ Cc il b jl 5 c i 2 b j2 9 


G. b. 1 = <c.b.> 
m~jn J vj 


and 



<c..d> 


where to n and w are taken as impulses. 
Pi 9 


Table 2 Natural Frequencies and Damping Ratios of Eight Cases 
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Figure 1 Tracking Performance Versus Cases 
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Figure 2 Cooper-Harper Ratings 
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The remaining piece of information about these configurations relates 
to the magnitude of the participation of the modes in the attitude response 
of the vehicle. This is measured by the relative magnitudes of the resi- 
dues of the modes. 

Now four vehicle responses known to be significant in longitudinal 
dynamics are flight path angle y, vertical acceleration at the cockpit n z , 
rigid-body pitch attitude o, and pitch rate 6, In addition, the local 
indicated pitch attitude and pitch rate (Op6j) are also important for flexibl 
vehicles. Local indicated attitude differs from rigid attitude in that 
the former includes the local elastic slope of the structure relative to 
the rigid-body axis. This local attitude angle is referred to as "indica- 
ted" since this is the quantity measured by a sensor placed at the 
fuselage station in question* 

Shown in Figs. 3-6 are the residues fur the two rigid-body modes 
(Phugoid and short-period) and as well as four aeroelastic vehicle modes 
(E1-E4) for four of Yen's eight configurations. Two of thesemodes (El 
and E3) are the two modes included in the configurations simulated by 
Yen (who did not include the two less important modes E2 and E4 in his 
math model). All these modes are identified from their eigenvectors. 

Now the problem with configuration 3 is clear. Unlike the other three 
configurations, the rigid-body attitude response, 0, is dominated by 
elastic mode El rather than by the rigid-body short-period model This is 
evident from the relative magnitudes of the modal residues in the theta 
responses. Note here that the input is an impulse passed through the 
pilot's lag model, or the actual input if the pilot attempted to input 
an impulse. From this result, it is clear that any attempts by the pilot 
to track a commanded attitude would be amost impossible. Also, character- 
izing the configuration handling qualities in terms of the rigid-body 
eigenvalues is meaningless. 
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(DUE 10 PILOT INPUTS) 
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Included as Appendix B are the residues for all eight configurations 
for all six responses to pilot-filtered impulses, Also giver, in Appendix 
B are the residues for the same vehicle responses to atmospheric turbulence 
input, This input consists of an impulse passed through a Dryden gust model, 
and thereby provides indications of modal contributions to the vehicle's gust 
response. Upon comparing the "gust response residues" to the "pilot response 
residues" we note for example that the situation regarding o residues on con- 
figuration 3 and 4 H very different. In the case of "gust response" the 
attitude o response is dominated by mode P.3 for configuration 4, while that 
was not the case for the "pilot response" for this configuration. This 
would indicate that a task involving attitude response in turbulence would 
be very difficult for configuration 4, while this configuration was much 
better than configuration 3 in attitude tracking, which doesn't involve 
vehicle gust response, finally note the large participation of the elastic 
modes in the acceleration response (n z ) of all configurations for both types 
of inputs. * 

With regard to the lateral -directional dynamics, two pilot Inputs, 
aileron and rudder, are involved of course, and the significant vehicle 
response variables must be selected before residues may be found. Impor- 
tant aircraft responses in this axis include sideslip (b) and bank 4) 
angles, roll (p) and yaw (r) rates, and lateral acceleration at the cock- 
pit ly In addition to these, experimental evidence points to the fact 
that the rate of onset of acceleration may be more Important than the mag- 
nitude of lateral acceleration in ride and flying qualities. Therefore, 
an additional response is included in the analysis, and it is defined as 
"lateral jerk," or diydt. Finally, both the rigid-body portion of these 
response variables as well as the total, or indicated responses, which 
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include rigid plus elastic contributions are included in the response 
variables considered. The "total" responses are again referred to as 
indicated since these variables would be measured by the appropriate 
sensor, and include total (rigid plus elastic) response. Consequently, 
the vector of responses considered in the lateral directional analysis 
is taken as 


V I.D. 3 O r ig’ Prig 1 r rig> n Y Hg » n Y Hg * *1' p r r r n Yj’ n Y^ 

The vector of control inputs is u T « $ R ]. And finally, pilot lags 

on rudder and aileron are included in the dynamiv ,iodel , along with a 
Dryden gust model appropriate for this axis. 

For the baseline B-l type vehicle, the eigenvalues are listed below* 
while the residues for pilot and gust responses are shown in Figs. 7-9. 
Figs. 7 and 8 include the response residues for an impulse input through 
the pilot (lag) aileron and rudder, respectively. Fig. 9 includes the 
gust response residues. The significance of the elastic modes in many of 
these aircraft responses is quite evident, even for the baseline configu- 
ration. 

Attention is now turning to the parametric variation of in-vacuo 
frequencies for the lateral -directional model, and thereby generate a 
family of "more flexible" configurations, as in the longitudinal case. 
Analysis will proceed then on these configurations, along with planning 
a fixed-base simulation. 


MAGNITUDE MAGNITUDE MAGNITUDE 


ORIGINAL PAGE 19 
OF POOR QUALITY 


16 


RESIDUE MAGNITUDES 

(DUE TO PILOT INPUTS) 

CASE 1 - "DA 


RESIDUES FOR ERCH OUTPUT ARE NORMALIZED SO 
THAT THEIR SUM IS 1.0 


PHI PHI— I 



* 

a 1 - 0 -] 
g B _ 






£ .0 - 

CD 

CE n 

51 ,0 - 1 





RL 

DR 

SP El E2 E3 EH 
MODES 

RL 

DR 

■■» ■ ■ . ... 

SP El E2 E3 EH 
MODES 


PHI DOT 


PHI DOT- 1 


1.0 

.5 


.0 


I 1 


RL DR SP El E2 E3 E4 
MODES 


ul.O 


ZD 

I— 


Z 

CD 

a: 

sz 


.5 


£ .0 


RL DR SP El E2 E3 EH 
MODES 


1 .0 

.5 


.0 


YAW RATE 


RL DR SP El E2 E3 EH 
MODES 


uul.O 

ZD 

t .5 

z 

CD 

CH n 

.0 


YAW RATE- I 


RL DR SP El E2 E3 EH 
MODES 


Figure 7 






MAGNITUDE MAGNITUDE 


17 


r WGINAL PAGE IS 
OF POOR QUALITY 


NY 



NY— I 


£ 1*° 
Z3 


Z 

CD 

cn 


.5 

.0 


RL DR SP El E2 E3 E4 
MODES 


RL DR SP El ES E3 E4 
MODES 


NY DOT 


NY DOT- I 



gl.O-n 

S B _ 




n 

L U L_ 

■z. 

CD 

£ .CH 



n 




RL DR SP El E2 E3 E4 RL DR SP El E2 E3 E4 

MODES MODES 


Fig' re 7 , cont'd. 





MAGNITUDE MAGNITUDE MAGNITUDE 


ORIGINAL PAGE IS 
OF POOR QUALITY 


18 


RESIDUE MAGNITUDES 


(DUE TO PILOT INPUTS) 

CASE 1 - DR 


RESIDUES FOR EACH OUTPUT ARE NORMALIZED SO 
THAT THEIR SUM IS 1.0 


PHI 


PHI — I 



RL DR SP El E2 E3 E4 
MODES 


uj 1 .0 

C3 


SC 

CO 

«r 

s: 


.5 

> q . . 




RL DR SP El E2 E3 E4 
MODES 


PHI DOT 


PHI DOT- 1 


. - — , 




n=L 

a - 


RL DR SP El E2 E3 E4 
MODES 


w 1.0 


2 

ca 

cn 

xz 


.5 4 

.0 


RL DR SP El E2 E3 E4 
MODES 


YRW RRTE 


YRU RRTE- I 


1 . 0 ' 

.5 

.0 


RL DR SP El E2 E3 E4 
MODES 


uj 1.0 


1 r* 

r-4 .b 

Z 

CD 

cn 

^ .0 


RL DR SP El E2 E3 E4 
MODES 


Figure 8 







MAGNITUDE MAGNITUDE 


19 


ORIGINAL PAGE IS 
Or POOR QUALITY 


NY 


1.0 

.5 

.0 


~i 


izu 


RL DR SP El Eg E3 EH 
MOOES 


NY— I 


iu 1 .0 

o 


H- 

'JC 

CD 

tr 


.6 

.0 





tltlWOTCM 




RL OR SP El Eg E3 EH 
MODES 


NY DOT NY DOT— I 





u.j 1.0-, 

Cl 

p 




r _ m 


r-n — 

£ .5“ 

CD 

cc 

s: ,0-i 

__ r~i 




RL DR SP El E2 E3 EH ' RL DR SP El Eg E3 EH 

MODES MODES 


Figure 8, 


cont'd 





MAGNITUDE MAGNITUDE MAGNITUDE 


20 


ORIGINAL PAGE IS 
OF POOR QUALITY 

RESIDUE MRGNITUDES 

(DUE TO DISTURBRNCE) 

CASE 1 - GUST 


tm KVig ?T T " RE N ™ RL,ZE ° 


SO 



RL DR SR El E2 E3 EH 
MODES 


PH I — I 


UJ 1 .0 

D 

Z) 

S -5 

CD 

01 r, 

= .0 


RL DR SP El £2 E3 EH 
MODES 



PHI D0T-I 


u 1.0 

~ .5 

z 

CD 

CE _ 

= .0 


RL DR SP El E2 E3 EH 
MODES 


RL DR SP El E2 E3 EH 
MODES 


YflW RRTE 


YflW RFITE-I 


1.0 


.5 


.0 

RL DR SP El E2 E3 EH 
MODES 



ul.O 


CD 

CE 


.5 

.0 


RL DR SP El E2 E3 EH 
MODES 


Figure 9 








ORIGINAL PAGE IS 
OF POOR QUALITY 






Table III 

Lateral -Directional Eigenvalues (Baseline) 


Spiral, - .050 Elastic 1 , -1 .959 ± 16.389j 

Roll, -2.191 Elastic 2, -0.570 ± 15.599J 

Dutch roll, -0. 159+2. 002j Elastic 3, -0.286 + 20.288j 

Elastic 4, -3.753 ± 23.849j 


3. Personnel and Plans 

Three graduate students are currently supported on tMs grant project. 
They are Mr. Marty Wasak, an M.S„ student currently responsible for the 
dynamic analysis of the longitudinal axis. Mr. Wasak will also perform 
a Neal -Smith analysis of these configurations via an optimal control pilot 
model, as developed in research projects here at Purdue. This work will 
constitute Mr. Wasak's M.S. thesis to be completed in Aug. 1983. He has 
just indicated that he plans to remain after August to pursue doctoral 
studies, and will therefore continue to be supported by this grant. 

Mr. Frank Leban is the second M.S. student supported by the grant, 
and is currently completing the computer software for the lateral -directional 
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dynamic analysis. This software should be completed by May, 1983, and 
results from this analysis will thereby be available for the next status 
report, 

The third student recently starting on the project is Mr. John 
Davidson, also an M,S. student. Mr. Davidson began investigating a modal 
control technique considered applicable to structural mode control. This 
control synthesis method is based on a method described by Moore, and 
extended by Cunningham of Honneywell. As yet, it had not been applied 
to flexible vehicle control, but it appears to be a prime candidate for 
shaping the response of such vehicles. In application of this method, 
optimum handling characteristics, will be the goal. 

A fourth student, Mr. William Smith, will also contribute to this 
research project, though as an Air Force officer he will not be financially 
supported by the grant. Mr. Smith's responsibilities are expected to 
focus on the fixed-base simulation of the lateral -directional dynamics of 

t 

a family of flexible vehicle configurations analogous to the eight confi- 
gurations on which the longitudinal analysis has been performed. It is 
hoped that these simulations will begin this summer. Roll tracking will 
be the task investigated, and it is to be noted that the Purdue's simula- 
tions will all be fixed base. 
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A Modal A ria 1 us is of Flexibl e Aircraft 
P unamics with Hand I inq Qual it ies I mol i cat ions 

by 


David K. Schmidt 

School of Aeronautics and Astronautics 
Purdue University 
West Lafayette, IN 47907 

Extended Abstract 


Int roduct ion 

The "flying qualities" or "handling qualities" or a 
flight vehicle refer to the ease or difficulty with which 
the pilot is able to control the vehicle in performing the 
various missions for which the system was designed. In 
this definition we see three important factors. They 
i tic lude 

The manual c ontrol lab .i 1 itu of a d unam ic system 

The use of sub j ect ive measures 

("ease. , .difficulties"^ 

The mult iol ic itu of missions and tasks . 

Inherent in the first two is the all important fact 
that a human is involved, hence psychological and physio- 
logical factors strongly come into play. Through extensive 
training of pilots, the attempt is made to minimize the 
deviation among these var iables and to minimize inherent 
human limitations. However, biases of opinion, likes, 
dislikes, and so forth, are still present. This fact, 
coupled with the ur>e of subjective measures make the 
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evaluation of flying qualities, not to mention the des i on 
problem, quite challenging, However, with a careful and 
systematic approach, both can be accomplished. 

Historical ly, the subjective opinions from highly 
trained test pilots on the performance and workload 
involved in performing the missions were correlated with 
the damping and frequency of the vehicle's rigid body 
modes. These parameters < damping and frequency) were then 
used to infer the handling qualities of future vehicles, 
or for handling qualities specif ic at ions . El 3 

This approach was very successful over the years, but 
in a very real way caused many of the problems < primarily 
conceptual) experienced recently with modern, more complex 
vehicles. This damping and frequency became sunonumous 
with the vehicle's handling qualities. That is, they 
became thought of as the handl ino q ual it ies rather than 
two parameters that affect the handl ino qualities ! 

Here is the key point - the handling qualities of a 
vehicle are not determined by the damping and frequency in 
a transfer function, but by the dynamic behavior, the time 
res ponse , of the vehicle system. This simple fact 
explains why the classic approach was successful for 
rigid, unaugmented aircraft, while high— order dynamic 
effects < aeroelast ic or control system dynamics) have been 
found to so significantly affect the pilot opinion as to 
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render the .approach useless. 

The purpose of the work to be presented in this paper 
is to address the quest ion, "Can dynamic aeroelastic 
effects significantly affect aircraft handling qualities,, 
and if so,, how?" By answering this question we are then 
able to consider the enhancement of the handling quali- 
ties, if required, through active control or structural 
design. 


To move toward answering this question, we will dis- 
cuss an analytical approach for evaluating a subject 
vehicle*' s dynamics, and apply the method to a spectrum of 
generic vehicles for which experimental data has been 
obtained. The complete methodology includes not only open 
loop analyses of the aircraft dynamics, but also an 
ana 1 y t i c a 1 pilot /veh ic 1 e ana 1 ys is to be present ed . 

M odal, hns hjsi.s Conce pts 

In the transfer function for the shorf -period approx- 
imation for longitudinal aircraft dynamics, the numerator 
depends on the same dimensional stability derivatives as 
those in the characteristic equation. CSU So the transfer 
function does not involve four independent parameters <two 
poles, a zero, and gain> but rather more like three. 
Therefore, once the stick gain is fixed, and short period 
damping and frequency are specified, very little 
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variability remains, As 49 result, not# that th# short 
Period time response i* essent ial ly specified .ns well 
^assuming the JiMar model to b# valid>, The point is 
that in this case, specifying the rigid body up^aJL charac*-* 
t eristics nearly completely specifies, th& as 

we 1 1 , 


Clearly,, for systems for which a simple, rigid-body 
transfer function is a poor model, the parameters in that 
transfer function cannot guarantee good or bad handling 
qualities, for they cannot alone determine the system's 
time response . This is? our first important point - p arjam- 
pterj£ in a vehicle dynnmix model 1£ May 
piqvnfic^ntly affect th& r^f&Mfi'SE tin\& re^ponffja. Note 
that this is an operational definition for a systems dom- 
inant modes, but is more general . 

Let us now turn to horn parameters in a dynamical 
model affect time response, Again, assuming a linear 
model to be valid (or valid enough to gain insight >, the 
impulse response may be found from the partial fraction 
expansion of the transfer function. Taking the important- 
pitch attitude response, for example. 


0<S> 


n 

E 

i»l 



u»here 


5 

n « order of the system 
« ith residue X^ 

X^ a ith system pole or eigenvalue 

Then the time response is 

# n -X . t 

8<t> * s R.e 1 
A»1 

Now for a classical rigid vehicle a good model mas 
obtained with n simply equal to two <3>, or a short period 
approximation. Now n may have to be larger than two for a 
good model , 
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Here rather than the time response expressed in terms 
of eigenvalues., it will convenient to the use the system's 
modes, So we may express the same time response as 

. N -o.t 

0<t.'> = S IR. le cos Co.t + T , > 

i=i 1 1 1 


where N *= number of system modes 


± jeix - X. X,. < conjugate pair of eigenvalues* 


T i J 

R^ = IR^I e < residue associated with root X^ > 


* 
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So the time response is completely determined by the 

real and imaginary parts of the eigenvalues <0^ and or 

bu £ . and o > . And the magnitude and phase of the resi- 
i 

dues <IR. I and •?.> associated with each mode. Or not only 
the modal damping and frequencies are involved in defining 
the time response# but the residues as well. Furthermore# 
in general these are the minimum number of var iab les 
required — AN. 


Returning for a moment to the classic# rigid body 
case# the use of only the short period damping and fre- 
quency in specifying the handling qualities is justifiable 
since the residues depend on the interaction between the 
poles and zeros of the transfer function# and as cited 
previously# the numerator character is tics are not indepen- 
dent of the denominator in a rigid aircraft transfer func- 
tion. Therefore# the variation in residues was limited. 
Note also# that the rigid-body assumption implies the 
r esidues for the infinity of elastic modes actually 
present in all vehicles are small compared to those of the 
rigid-body modes. 


This is certainly not so with the introduction of 
significant aeroelastic effects and control system dynam- 
ics# for example. The modal characteristics are still 
fundamental in their affect on time response# but the 
"numerator dynamics" or the interaction between the poles 
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and zeros are significant in terms of their effect on the 
res i dues . The zero's locations in themselves are not 
important , but rather how the position affects the resi- 
dues , 


Let's consider a simple example involving the short - 
p»eriod plus a single elastic mode,, for which the transfer 
furic t i on iti i gh t be 

e<s> 

A<s> ' 2 22 2 

* <s rf + 2c> <0 s + Xs* + 2£„ o, s + «?> 

sp sp sp '11 1 

be 

<S l t + T 1 1 
or 


The impulse (time) response of this vehicle may 
written as 


~<j t 


I- Uj -n u> V * T -r irs,. 

i*p sp sp 1 




9<t> — IR _le “' r ' cos<5 t + y > + IR, le co; 


K<s+z> <s^-+as+b> 


with 


e<t> 


0 _, . ■ , + 6 ,- , 

Rigid Flex 




co 


1 


<0 


1 
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Now, the affect of the f 1 ex i b 1 e mode on the time response 
consists of three components. 


The coupling affect on the short period modal 

characteristics < changing »7 _ and •» > 

sp sp 

The coupling affect on the short period residues 


<c hanging IR I and Y__> 


The time response of the mode itself <8,-, > 

Flex 


This leads to a second important point . The effects 
of higher-order modes mau be quart if ied in t erms of these 
effects on the dominant < u sual lu rioid bodu > modal charac- 
teristics and residues , as wel 1 as the time response of 
the h i gher- modes t hemse 1 ves . 


The above presentation u»as a review of known facts, 
but demonstrates concepts upon which we have developed our 
ana lytical approach . For more realistic dynamics <i.e. 
higher-order> , numerical analysis using mult i -variable 
systems techniques is appropriate. We assume a linear 
state-variable represent at ion for the vehicle is avail- 
able, and although the mult i- input , multi -out put case 
<e.g. 8 and n_> will be treated, consider for now the case 
above, or 8/ . We have 


x = Ax + B£_ 


6<t> = Ox 
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Now diagonalizing the system via the model matrix T con- 
sisting of the system eigenvectors V^, or 


T = CV.V 0 . . , V 3 

Id! n 


one obtains 


q — Aq + T * B.Sj- 

the vehicle dynamics expressed in modal coordinates C33. 
With the pitch rate expressed in terms of the system modes 
as 


8 < t > = CTq 


where A = diagonal matrix of eigenvalues., X^ = ± jo; 

I* * 

— 1 I 1 1 

T B = Modal controllability matrix, .• , 

b 

I n I 

CT — Modal observability matrix, . . .£ , 

” 1 TlJ 

And the partial fraction expansion of the system's impulse 
response is simply 


n 


e i B i 


— & . 

. , s-X . 

x=l 1 


e<s> 
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obtained directly from Laplace t ransf orm i ng the above 
equations. So me have the residues simply as the pro- 
duct £\ , available by inspection from the sustain s modal 
re present .at i on . The three effects of a higher-order mode 
discussed previously may therefore be found numerically. 

The ftna 1 us i s Approach - Overv i ew 

Although all the above ideas are available from 
linear systems theory, blind use of these ideas mill lead 
to incorrect answers to our fundamental questions cited at 
the outset . The straightforward modal analysis of this 
linearised math model of the aircraft dynamics alone 
ignores the salient characterise ics of the actual distur- 
bances exciting these dynamics in actual flight condi- 
tions. Two such disturbances of interest here are the 
pilot's control input and atmosphere turbulence. Neither 
of these may be treated as "white” random processes, for 
example. The significant characterist ics of these "excit- 
ing subsystems" must be taken into account. How this is 
accomplished will thoroughly explained in the complete 
paper . 

The counterpart to this "open- loop" analysis is the 
p i 1 ot /veh icle analysis, performed us i ng an opt i ma 1 - 
control— theoret ic model of the pilot's control loop clo- 
sures. We will show, unlike the reduced-order pilot model 
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of Swaim C43, how our full-order pilot model . consistent 
with that of Ref. CS!, may be used to obtain good correla- 
tion between tracking error statistic* obtained analyti- 
cally and those obtained from fixed base simulation. 
Thorough this pilot modeling approach we are able to infer 
a different piloting strategy rather than a different 
pilot "internal model" C43 in the control of flexible air- 
craft dynamics, and this supported by both simulation 
results and the open- loop aircraft analysis presented pre- 
viously . 


Sample Results 

Shown in Figure X are the rms attitude tracking 
errors obtained from previous fixed base simulation of 
several generic flexible vehicles C63 . The onlu difference 
between these vehicles is their structural stiffness, 
quantified in terms of the in- vacuo vibration modal fre- 
quencies. The effects of flexibility in the vehicle 
dynamics are clearly evident in the increased tracking 
error, particularly in configuration C^. What is not 
clear is why configuration C^, an equally flexible confi- 
guration to in terms of elastic modal frequencies, 
resulted in significantly better tracking. The modal 
analysis technique out 1 ined above and to be described com- 
pletely in the paper will reveal the answer to this para- 
dox. In addition, how the full-order pilot model is 


12 


ORIGINAL PAGE IS 
OF POOR QUALITY 


applied to achieve the correlation between analytical and 
experimental results shown in the figure mill also be dis- 
cussed. These are just four of eight vehicle configura- 
tions that mill be presented/ with analytical results com- 
pared to experimental . In all cases, excellent .correla- 
tion is obtained. Finally, the modal analysis method is 
used to demonstrate the affects of flexibility on a 
variety of vehicle responses, in addition to pitch atti- 
tude, and these results mill be presented as well. 
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